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^ : 1 Introduction 

Q ■ EROS has discovered ~ 70 SNs during 8 periods partially dedicated to a 

SN search. The average discovery rate is 1 SN per 2 hours of observation, 
^ I corresponding to expectation. Before 1999, the SNs that we discovered were 

publicly available on the web and reported in lAU telegrams. The low sci- 
entific return of this procedure convinced the EROS collaboration to work 
within a larger collaboration. Nevertheless, EROS has taken advantage of the 
homogeneity of its SNIa sample to obtain a SNIa explosion rate at z ~ 0.1|jl|. 

EROS has participated in the large campaign of February-March 1999 
lead by the SCP, which aimed at measuring ~ 10 SNIa's discovered before 
or at maximum, including photometric and spectroscopic follow-up. EROS 
contributed an homogeneous sample of 8 SNIa's discovered before maximum 
(of which 7 were monitored). Our procedure was as follows: The first SN 
detection was made on the day following the night during which 20 pictures 
(one square degree each) were taken. A confirmation picture was taken the 
following night, and the SN coordinates were then given to the coordinator. 
The complete process used to take less than 48 hours to provide a confirmed 
SN. This delay was clearly acceptable, if one considers the > 50% proportion 
of pre-maximum SNIa's delivered by EROS. A preliminary analysis of the 
photometric data collected for these SNIa's is described in detail in Nicolas 
Regnault's thesispl. 
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The EROS-2 microlensing search will end at the end of 2002. In this doc- 
ument, we investigate a new way of using the EROS telescope (the MARLY) 
after this date. 

2 The MARLY as a SN-photometer 

2.1 Today's status 

The MARLY present optics has a wide aperture, and the sampling is rela- 
tively poor (EROS was designed to maximise the number of monitored stars 
per pixel). The pixel size is 0.6 arcsec. A test, done with the heat sources 
of the dome switched off, showed that the dome contribution to the exper- 
imental seeing was 0.8 arcsec on the optical axis. As there is a significant 
degradation of the seeing off-axis and because of the effect of the heat sources, 
the median instrumental seeing is only 1.5 arcsec. 

2.2 Telescope configuration changes 

The Marly telescope had originally a F/8 aperture. For the EROS pro- 
gramme, a focal reducer has increased the aperture to F/5. With the original 
aperture, a simple setup with a one-CCD camera and a BVRI filter wheel 
could be installed on the telescope. 

We will also propose a more sophisticated setup to perform U band photom- 
etry in parallel with the other bands. 

2.3 Optical photometer performances 

In this section, we estimate the total time needed to fully monitor a SNIa 
with the MARLY telescope every other night in BVRI during 80 days after 
its discovery. 

2.3.1 Expression of the exposure time 

For a given passband w, the exposure time Tg^p needed to reach a relative 
photometric precision 5 on a pointlike object with apparent magnitude M^, 
depends on 

• the photoelectron flux $^(Mu, = 0) (in je'/s/m^) associated with a 
source of magnitude M^, = 0, taking into account the detector effi- 
ciency, 

• the magnitude of the background surface brightness fj^^ per arcsec^, 
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• the seeing. 



For the MARLY collecting surface, the calculation and the validity conditions 
of the following formula are given in annex A: 
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We will use this expression to perform a realistic simulation, taking into 
account real distributions for each parameter. 



2.3.2 Photoelectron signal expected from the SNIa's 

By definition, $^„(M^ = 0) depends on the spectrum of the measured point- 
like object. For SNIa's, we estimate the number of photoelectrons produced 
in the detector in each of the UBVRI passbands by using the following 
ingredients: 

• The photon spectra and the magnitudes of two typical SNIa's: 
SN1990N, 7 days before and 7 days after maximum^], and SN1994D, 
3 days after maximum 

• The atmospheric transmission curve taken at airmass=l. 

• The Bessel UBVRI theoretical transmission curves. 

• The LBNL red-enhanced CCD spectral response. 

• The transmission of the MARLY F/8 optics, which has been assumed to 
be uniform in wavelength, equal to 0.87 per reflector, i.e. 0.87^ = 0.76. 

Fig. ID shows the various transmission curves and the SN1994D spectrum. 
From this input, we find that a SNIa with Mb — Mv = 0, and a spectrum 
typical of the one emitted within +/-7 days from maximum, should produce 
a flux of photoelectrons in the detector given by the upper part of Table 
1. Taking into account the measured colour indices of SN1994D, we get 
a correspondance between the zero Myj magnitudes and the photoelectron 
fluxes in UBVRI (lower part of Table 1) []. 

"'^ These estimates are compatible with the observed signal on SNIa images taken with 
the Danish telescope. But only a rough check could be done because of differences in the 
filter transmissions and CCD efficiency. 
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Figure 1: Upper panel: transmission of the UBVRI filters (coloured lines), 
the atmosphere (black thick), the optics (black standard), and quantum effi- 
ciency of the CCD (orange thick) as a function of the wavelength. 
Lower panel: measured SN1994D spectrum in photons/rn^ / s/ A, 3 days af- 
ter maximum (black thick line), residual spectras after atmospheric, instru- 
mental and filter transmissions, and conversion in the CCD in jm} j sj A 
(coloured lines). 
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Table 1: The expected total 76 fluxes produced by the LBNL red-enhanced 
CCD at Airmass=l for a SNIa at maximum in B. The collecting area of the 
MARLY IS S = 0.57 m\ 
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Table 2: Sky brightness in magnitudes for different passbands, as a function 
of the number of days from the New Moon; range and average of the host 
galaxy brightness. All magnitudes are the ones at the top of the atmosphere. 
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2.4 SN photometry simulation 

We have computed the total time spent in the photometric follow-up using a 
simple SN light curve simulation. Wc suppose all SNe to lie at z = 0.05. The 
whole calculation is done for airmass=l. The SN discovery age is supposed 
to range between 5 and 10 days before the time of B maximum with a flat 
distribution. The lunar phase is also randomly chosen. 

2.4.1 SN environnement - Expected photoelectron background 

The background surface magnitude due to both the sky brightness and the 
host galaxy surface brightness in the w passband is given by: 

When the background flux is negligible (i.e. is large, or if the SN is 
bright), the last factor in the expression of Texp is 1, and the formula can be 
simplified. 

• Night sky brightness ^^'^^ 

The sky brightness depends on the Moon phase and also on the atmo- 
spheric conditions. We adopt the average values used for the prepara- 
tion of the ESO proposals (Table 2, upper part). 
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Host Surface Brightness |J,g(mag/arcsec ) 

Figure 2: Distribution of the host galaxy 'B-magnitude in the EROS SNIa 
discovery sample. 



• Host galaxy brightness /i^ 

The background hght from the host galaxy depends on the position of 
the SN. If the SN is in the disk, far from the bulge, then the average 
surface brightness of the galaxy in B is fi^'' = 22 A 0. If the SN is 
located in the bulge, the average is /Xg"' ~ 20. The SNIa's referred to in 
the EROS paper |I| have a local host surface brightness ranging between 
^gii _ -j^g ^gi _ QT^(^ more than 80% have > 20 (see Fig. 

0). Following ^ we adopt < /i^' — /ig"' >~ and < — /xf^' >~ 1 
as typical values, and we extrapolate for V and I. In the following, we 
will assume that the SNIa's that we will follow-up will belong to such 
a sample. 

2.4.2 SN light curve simulation 

Light curves of SNIa's at z = 0.05 are generated according to the observed 
frequency of various templates described in 0. 

SNIa light curves show some diversity. We use B,V,R and I templates 
obtained in [0] to describe the different luminosity evolutions with time. We 
used those templates from reference [@ for which early data is available 0. 

2SN 91T, 91bg, 92A, 94D, 95D, 96X, 98bu, 98de 
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Figure 3: SN Light curves (in mag) obtained in the simulation. 
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Figure 4: Distribution of the seeing at La Silla (in arcsec). 



The absolute peak magnitude attributed to each template is based on the 
standardisation relation determined in [Q]. 

The program randomly generates a template according to the frequency 
distribution which came out of in the best fit analysis of 0. The SN max- 
imal luminosity is smeared using a normal law with 0.1 mag dispersion (in 
each band). For each SN, a host surface brightness is generated according 
to the distribution shown in Fig. ||. Once a discovery age, lunar phase, host 
brightness and SN subtype have been determined, the program follows the 
SN luminosity evolution for 70 days after B maximum. Photometric measure- 
ments are assumed to take place every other night, and a 70% good weather 
is assumed. In a bad weather case, the new measurement is attempted again 
the next day. 1000 light curves obtained this way are shown in figure |^. 
We associate a random seeing to each measurement, using the true seeing 
distribution of La Silla (see Fig. |); an instrumental contribution of 0.8", as 
measured on the MARLY optical axis, is quadratically added to this seeing. 
The exposure time necessary to reach a 2% photometric precision is then 
computed for each measurement. These times are shown in figure |^. These 
exposure calculations are done for the MARLY telescope. The scaling to be 
performed in order to get results for another telescope (at the same site) or 
another resolution is straightforward. 
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Figure 5: Logio{Texp) needed to reach 2% photometric precision versus the 
SN age (in days). Exposures taken around full moon (^±4 days) are shown in 
red. The horizontal lines illustrates the 1 hour cutoff used in this study. 
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Figure 6: Logio{Texp) needed to reach 2% photometric precision versus the 
Moon phase (in days). The red dots correspond to SN age > 15 days. For 
these late measurements, the exposure times are much longer and strongly 
depend on the Moon phase. The horizontal lines illustrates the 1 hour cutoff 
used in this study. 
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Figure 7: Distribution of the total follow-up time per simulated SNIa, in 
hours. Each entry corresponds to the sum of the exposure times of about 40 
measurements with 2% photometric precision. 
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2.5 The "cost" of a SNIa 

Fig. 1^ and |] give the distribution of exposure times as a function of SN age 
and Moon phase respectively. Long exposure times occur when the SNIa 
fades, typically from 15 days after maximum (red dots in Fig. |^). The 
exposure then strongly depends on the sky brightness in B and V. 

We estimate the telescope time needed per SNIa by adding the exposure 
times of all the measurements (about 40 measurements per SNIa and per 
passband). An overhead of 60 s is added to each exposure. For long expo- 
sures, which will be fractionned, we further add 60 s of overhead to every 
sequence of 20 minutes of observation. We choose to truncate the exposures 
at 3600 seconds. Such a truncation affects 7% of the measurements and 
takes place only in the late part of the SN light curves {age > 30 days). 
The photometric resolution is then slightly downgraded (from 2% to 4% in 
the worst case). Fig. |^ shows the histogram of the total telescope time for 
each passband, for our sample of 1000 simulated SNIa's. The outcome of 
this study is that the follow-up of a SNIa, involving about 40 measurements 
regularly spaced within a 80 day period, needs an average telescope time of: 

• 9.5 hours in B, 

• 4.5 hours in V, 

• 3.8 hours in R, 

• and 5.3 hours in I, 

Table 3 gives the results of the same calculations done for other CCD effi- 
ciencies (see Fig. H). 

Table 3: Total follow-up time (hours) for 4 CCD efficiencies, in the BVRI 
bands 



CCD 


B 


V 


R 


I 


Total 


Berkeley 


9.5h 


4.5h 


3.8h 


5.3h 


23.1 


SITe 


8.6h 


4.2h 


3.7h 


6.2h 


22.7 


EEV 


10.6h 


4.6h 


4.2h 


7. Oh 


26.4 


Kodak 


14.1h 


6.1h 


4.9h 


9.4h 


34.5 
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Figure 8: CCD efficiencies used in the calculation. 



2.6 Strategy 

In order to minimize the cost of the measurements taken during the late 
section of a SNIa (when faintest) we could consider some improvement to 
the very simple strategy described above: 

• As long as the SNIa (at z = 0.05) is close to its maximum (~ 17.5), one 
can neglect the sky brightness in BVRI and the host galaxy luminosity 
in a first approximation. This is no longer the case when the SN fades. 
Then, when the Moon is bright, we should favour data taking on the 
brightest SNIa's, and devote the dark fraction of each night to the 
observation of the faintest ones. 

• As the luminosity variations of "old" SNIa's are slow, relaxing the 
sampling rate should have less impact than around the maximum. 

• Lastly, we could downgrade the photometric precision when the esti- 
mated exposure time is long. The impact of such a choice should be 
studied. 
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Figure 9: Exposures needed to reach a 2% photometric precision in U versus 
the SN age (left) and Moon phase (right) with the Berkeley CCD. The colors 
are drawn as for figures ^ and ||. The horizontal lines correspond to 1 hour. 

2.7 An additional dedicated camera for U 

The U band is of special interest, because the detailed study of the light 
curves of the nearby SNIa's will be useful for the I band study of high redshift 
SNIa's {z > 1). 

For this passband, only two SNIa templates are available. The first order 
calculation performed in the first version of this note (May 2001) showed that 
more than 60 hours were needed to follow-up a SNIa in U band, with our 
simple strategy. Figure || shows the individual exposures in U, computed as 
for the other filters for the Berkeley CCD, as a function of the SN age and 
lunar phase. 

Since unrealistically long exposures are needed for the U band, we have 
considered the possibility to take images in parallel in the U band and in 
BVRI with a dichroic beam splitter. A possible setup to achieve this is 
displayed on figure [1^. Preliminary discussions with optics engineers from 
Observatoire de Marseille confirmed the technical feasibility of this setup. 

In this setup the CCD detectors can be separately optimized for the U 
and the BVRI bands. The data taking could proceed as follows. Consecutive 
exposures in BVRI are chosen so as to ensure a 2% photometric resolution. 
The exposure in U lasts until the end of this cycle. To evaluate the viability 
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Figure 10: A possible setup for a photometer integrating in U and in BVRI 
in parallel 



of this scheme, we have estimated the resulting photometric precision in U. 
We have assumed that the UV-enhanced SITe CCD will be used for this 
channel (see figure H). Figure |ll| shows these individual resolutions as a 
function of the SN age and lunar phase. 

The moon has clearly a strong degrading effect (near full moon) ; late 
epochs (after day 25) are also attributed a worse resolution. 

Figure |T2| shows the resolutions achieved with the naive strategy discussed 
here. Almost all measurements have a resolution better than ^ 10%. The 
average resolution is ~ 3.5%. A more subtle observing strategy could be 
sought to optimize the achieved resolutions versus the observing time. For 
example, one could decide to sometimes impose the BVRI exposures to 
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Figure 11: Photometric resolutions achieved when exposing in U with a SITe 
UV- enhanced device, in parallel with BVRI and the LBL CCD, as a function 
of SN age and Moon phase. Red points indicate measurements scheduled near 
full Moon (left) or at a SN older than 15 days (right). 

reach a fixed, better resolution in U. 

The numbers obtained here in the U passband are somewhat uncertain 
because of the large fluctuations of the atmospheric absorption. However, 
since exposures in U and other bands are simultaneous, color dependent 
absorption corrections due to variations in atmospheric conditions should 
be controlled. Monitoring these absorption corrections will be easier with a 
larger CCD field ; 2A; x 2/c chips are at least necessary 

Finally, it is worth to mention that a factor of 2 in the photo-electron 
flux in U could be gained by removing the U filter. The effective passband 
(limited by the atmospheric absorption cutoff and the dichroic) would then be 
slightly modulated by the atmospheric conditions. This solution, still under 
study, might be considered only if the absorption is proved to be precisely 
monitored. 



■^In the i^/8 configuration of the MARLY, and with 15/i pixels, this corresponds to 
13' X 13' on the sky 
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Figure 12: Photometric resolutions achieved when exposing in U with a 
SITe UV- enhanced device in parallel with BVRI and the LBL CCD. This 
corresponds to a strategy aiming at a 2% resolution in BVRI. 



2.8 Maximum number of SNs that can be monitored 

Data will be taken sequentially with the BVRI filters. The total telescope 
time per SNIa will thus be: 9.5+4.5+3.8+5.3=23.1 hours, spread over a 
period of 80 days. The average duration of the astronomical night at La 
Silla is 9 hours, and the average clear sky fraction is 70%. Then the average 
available time per period of 80 days is 80x9x0.7=504 hours. Up to 22 SNIa's 
can be completely measured during this period of 80 days. Therefore the 
rate of 100 SNIa's per year completely measured in the UBVRI passbands 
is the capacity of the MARLY with the setup outlined in figure |10|. This rate 
could be made larger with a more subtle strategy as indicated above with a 
limited loss in sampling or in photometric precision. By using the tables and 
figures provided in this note, the reader can estimate the performance that 
would correspond to another strategy. 

2.9 Miscellaneous 

One should keep in mind that we have considered SNIa's with no reddening in 
our calculations. In the case of reddening, the estimates in the blue passbands 
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may be significantly affected. 

For the U band, airmass is a critical issue since the atmospheric absorp- 
tion is a dominant factor^. The signal from the SN and the host galaxy would 
be lowered by a factor 1.22 if the airmass were 1.5 instead of 1, and it is clear 
that the full Moon would strongly disturb the observations in this passband. 

Due to the presence of a dichroic splitter, the effective B band could be 
slightly narrower than the standard one. 

From figure |^ one may see that this setup is also sensitive in the Z band 
which could be added to our observation scheme in addition with BVRI. 
Nearby SNIa have not yet been thoroughly studied in this band which could 
be helpful to disentangle the SN reddening by the host galaxy's dust and 
the intrinsic SN colour dispersion. Finally, let us recall that the proposed 
configuration could be optimized to other photometric systems (e.g. ugriz). 

2.10 Necessary upgrades and schedule 

To use the MARLY as a photometer, according to our preliminary design, 
the following technical changes must be done: 

• The focal reducer has to be removed, and the F/8 optics has to be re- 
installed. A dichroic beam splitter at a wavelength near 4000 A should 
be installed. Moreover, the guiding system has to be adapted. 

• two CCD cameras have to be built : 

— A single CCD camera with a high performance CCD for BVRI, 

— A single CCD camera with a UV-enhanced CCD for U. 

The EROS project will stop taking data at the end of december 2002. De- 
signing and building the instrument can be largely decoupled from MARLY 
running and thus could be performed in France. Installing and testing the 
system on site should not take longer than one month, as in the ER0S2 
case. Provided a more general framework is found for this project, a startup 
as early as spring 2003 could take place. 

2.11 Expected performance 

Around 100 SNIa's per year at 2; ~ 0.05 could be fully followed-up in BVRI 
with a photometric precision 6 = 2%, and in the same time in U with a 3.5% 

^AU = 0. 44 /airmass and AB = 0.23/airmass at airmass=l, for a SNIa spectrum 
and the CCD spectral efficiency adopted here. 
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photometric precision on average. The exposure time scales with If the 
overhead time is small compared to the exposure time (if S < 2%), then the 
number of SNla's that can be monitored is ~ 25 x per year. U 6 > 2%, 

then the overhead time is not negligible and we have to use the following 
expression 



to correctly estimate the corresponding telescope time and the rate of mon- 
itored SNla's. We thus find that around 300 SNla's per year at z ~ 0.05 
could be fully foUowed-up with a, 5 — 4% photometric precision in BVRI. 

3 Conclusion 

Other possible uses (as a SN discoverer or with an integral field spectrometer) 
of the 1 meter dedicated MARLY telescope, based at La Silla, have also been 
evaluated in a more complete studiy available on request. 

Refurbished in order to simultaneously perform BVRI and U photome- 
try, the MARLY could be used for the complete follow-up (40 measurements 
spanning 80 days) of at least 100 SNla's per year at 2; ~ 0.05, with a photo- 
metric precision of 2% in BVRI and 3.5% on average in U. If the required 
photometric precision is only 4%, then about 300 SNla's can be foUowed-up 
each year. 

A project dedicated to the nearby SN study should integrate dedicated 
discovery and follow-up means, the latter providing spectrographic and/or 
photometric capabilities. With the equipment proposed in this note, the 
MARLY telescope could nicely complement such a project. We will would 
be happy to collaborate on any such project. 
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Annex A 

Expression of the optimal photometric preci- 
sion accessible with aperture 
photometry 

Let 

• Texp be the exposure time (in s), 

• F'^'^ be the integrated photoelectron flux from the SN in the telescope 
in photo-electrons per second {je~/s), 

• b he the photoelectron flux per pixel originating in the atmospheric 
emission and the host galaxy (in je~ /s) 

• p he the pixel size. 

Assuming that the PSF is a gaussian of standard deviation a — seeing/2.36: 



PSF{r)^- l^e-rV2.^ 

and that the pixel size is small with respect to the seeing., the time integrated 
photoelectron signal in a disk of radius r is given by 

5(r)=F^^xTe,j,x(l-e-'-'/2-^). 

During the same exposure time, the number of photoelectrons due to the 
atmosphere and the host galaxy in the disk of radius r is given by 



B{r) = bx T^^p X TT 



Assuming that we measure separately < B{r) > with a negligible uncertainty, 
then the fluctuation of the sum S{r) + B{r)— < B{r) > is given by 



a{rf ^ S{r) + B{r) + al x tt 



where ae is the readout noise per pixel. The photometric resolution is then 

S{r)- 



given by the ratio S = ^fl. The resolution 6 has its smallest value when 



^ = 0. Expressing this condition leads to the calculation of the optimal 
value of r, and the expression of the best resolution follows: 

5=[F^^xTe,J-5 xD, 



21 



R = nolsa^peak signai) 

1 2 3 4 5 



Figure 13: The degradation factor D of the resolution as a function of the 
ratio R of the noise per pixel to the maximum signal per pixel. 



where D is a "degradation" factor of the resolution, depending only on the 
noise to signal ratio in the pixel where the signal is maximal: 

bT^xp + total noise/pixel 

R = 



jpSNj' /2-ji z ^ peak flux of the SN /pixel 

For long exposures, the electronic readout noise will always be negligible 
compared to the fluctuations of the background brightness. Then R can be 
written: 



R= 1.13 



seeing 
1 arcsec 



2 



M-u 

10^ 



where M is the magnitude of the SN and fi is the total background magnitude 
per arcsec'^. The D factor variation with R is given in Fig. ^ ; we easily 
check that for R=0 (no background), the resolution is simply given by the 
gaussian fluctuation of the number of signal photoelectrons, as expected. The 
following formula is an excellent approximation (within 2%) of D, at least 
up to R=10: 

D ~ 0.27 + 1.5Vi? + 0.25. 

For a given passband w G {UBVRI}, F^^ is related to the magnitude 
of the SN and to the collecting surface S of the telescope through the relation: 

F^^ = $^(M^ = 0) X 5 X io-*^-/2.5^ 
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where ^w{Mw — 0) is the photoelectron flux in the w passband per second 
per collecting surface, for a SNIa with magnitude M.^, = 0. Then the 
exposure time Tf,xp needed to reach a photometric precision 6 is given by the 
expression: 



L>2 



- exp 



S X $^(M^ = 0) 



X 



5^ 



0.27+ LS-i 0.25 + 1.13 



seeing 
1 arcsec 



-\ 2 



10- 



This exposure time scales with the inverse of the telescope collecting power 
(surface x transmission) and with the inverse of the square of the requested 
resolution. The scaling with the seeing is more complex, and depends on the 
contrast between the signal (from the SN) and the background flux. 

For the MARLY telescope S ^ ^ x {dl^i^ - (P^^J = 0.57 where 



''prim. 



0.98 m is the primary mirror diameter and dg 



0.49 m is the 



diameter of the obstructing secondary mirror. Then T^xp can be expressed 
as follows: 



T ~ 

exp — 



(175s) 
[5{in% 



£«^(M^ 0) 

109 



n -1 



X 



Mro-17.5 

X 10 2.5 



0.27+ I.59WO.22 + 



seeing 
1 arcsec 



10 2.S 



When the background flux is negligible (i.e. if is large, or if the SN is 
bright), the last factor in the expression of Texp is 1, and the formula can be 
simplified. 
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